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Abstract. Measurements of CMB anisotropy are ideal experiments for discovering the non-trivial global 
topology of the universe. To evaluate the CMB anisotropy in multiply-connected compact cosmological 
models, one needs to compute eigenmodes of the Laplace-Beltrami operator. We numerically obtain the 
eigenmodes on a compact 3-hyperbolic space cataloged as m003(— 2, 3) in SnapPea 1 using the direct bound- 
ary element method, which enables one to simulate the CMB in multiply-connected compact models with 
high precision. The angular power spectra C;'s (2 < I < 18) are calculated using computed eigenmodes for 
5.4 < k < 10 and Gaussian random approximation for the expansion coefficients for 10 < k < 20. Assuming 
that the initial power spectrum is the Harrison-Zeldovich spectrum, the computed C;'s are consistent with 
the COBE data for 0.1 < fl a < 0.6. In low Q models, the large-angular fluctuations can be produced at 
ON ■ periods after the last scattering as the curvature perturbations decay in the curvature dominant era. 

Oh 1 introduction 

6 

The Einstein equation does not specify the global topology of the universe; therefore, there is no a priori 
reason to believe that the space-like hypersurface of the universe is simply-connected. If the space-like hy- 
persurface is multiply-connected on the scale of the horizon or less, there is a possibility of discovering the 
multiply-connectedness by the future astronomical observations. 

In recent years, there has been a great interest in properties of CMB anisotropy in multiply-connected cos- 
' mological models [2-8] . Precise measurements of the CMB anisotropy by the future satellite missions such as 
MAP and PLANCK may enable us to find the fingerprint of the multiply-connectedness in the CMB. There- 
fore, it is very important for us to simulate the CMB anisotropy in multiply-connected FRW models. 

Constraints on the topological identification scales using the COBE data have been obtained for some flat 
models with no cosmological constant [2,5]. The large-angular temperature fluctuations discovered by the 
COBE constrain the possible number of the copies of the fundamental domain inside the last scattering surface 
to less than~8 for these multiply-connected models. 

The mode functions for flat models can be analytically obtained; therefore, the angular power spectra are 
obtained straightforwardly. On the other hand, no closed analytic expression of the eigenmodes is known 
for compact hyperbolic (CH) spaces. Therefore the analysis of the CMB anisotropy in CH models has been 
considered to be quite difficult. To overcome the difficulty, the author proposed a numerical approach called 
the direct boundary element method (DBEM) for computing eigenmodes of the Laplace-Beltrami operator [9]. 
14 eigenmodes have been computed for a "small" CH space m003(-2,3) with volume 0.98139 in the SnapPea 
catalog and it is numerically found that the expansion coefficients behave as if they are random Gaussian 
numbers. 

We briefly describe the DBEM and the properties of the expansion coefficients which are used for computing 
the angular power spectra for low fi CH cosmological models. We calculate the angular power spectra for 



SnapPea is a computer program by Jeff Weeks for creating and studying CH spaces [1]. 



m003(-2,3) using computed eigenmodes with small k, and an approximate method for eigenmodes with large 
k, which are compared with the COBE data. 



NUMERICAL COMPUTATION OF EIGENMODES 



The advantage of the DBEM is that it enables one to compute the eigenfunctions much precisely than other 
methods as it does not rely on the variational principle and it uses an analytical fundamental solution, namely 
the free Green's function. 

Let us first consider the Helmholtz equation on a compact connected and simply-connected M-dimensional 
domain il in a simply-connected M-dimensional Riemannian manifold M. with appropriate periodic boundary 
conditions on the boundary dfl, 

(V 2 + fc>(x) = 0, (1) 

where V 2 = V J Vi, (i = 1,2, • • -,M), and V» is the covariant derivative operator defined on A4. A square- 
integrable function u is the solution of the Helmholtz equation if and only if 

ft[u(x), v(x)] = ((V 2 + k 2 ) u(x), v(x)) - 0, (2) 

where v is an arbitrary square-integrable function called weighted function and ( ) is defined as 

(a, b) = [ abdV. (3) 
Jn 

Then we put tt(x) into the form 

M 

u = ^2 u j<t > 3> ( 4 ) 

3 = 1 

where (f>j's are linearly independent square-integrable functions. The best approximate solution can be obtained 
by minimizing the residue function 1Z for a fixed weighted function v(x) by changing the coefficients Uj. 
Choosing the fundamental solution Ge(x, y) as the weighted function, Eq.(2) is transformed into a boundary 
integral equation, 

c(y)«(y)+/ G E (x,y)p-n^dS- f ^|^n* u ^/gdS = 0, (5) 
Jdo. ax Jan ° x 

where c(y) = 1 for y € Vt and c(y) = ^ for y G dil. Since f2 is compact, the discrete eigenvalues are represented 
as fco = < ki < k2, ■ ■ .. The mode k = k\ is the most important one that has the longest " wavelength" defined 
as 2ir/k. The author succeeded in computing 14 eigenmodes for k < 10 on m003(-2,3) [9]. fci is numerically 
found to be 5.4. 

From now on, we limit our consideration to CH spaces whose universal covering space is 7i 3 . We set 
the c urvatur e radius of H 3 to 1 without loss of generality. For convenience, we expand the eigenmodes u v , 
(iy = Vk 2 — 1) on the CH space in terms of eigenmodes X„;Yj m 's on Ti 3 , 

u v = ^2 X v i {x)Yi m (6,(j)), (6) 

Ira 

where X u i is a radial eigenfunction, Yj m is a spherical harmonic and £„/ m is an expansion coefficient. It has 
been numerically found that ^; m 's (for I < 19 and v < 9.94) behave as if they are random Gaussian numbers 
for m003(-2,3) , which is consistent with the prediction by random matrix theory. Random Gaussian behavior 
is also confirmed in a 2-dimensional CH space [10]. Note that some properties of a quantum system whose 
classical counterpart is a chaotic system can be explained by random matrix theory [11]. It has also been 
numerically found that the variance of Culm's is proportional to v~ 2 . Using these properties, one can compute 
the approximate contribution from "highly-excited" modes (fc»l) to the angular power spectrum. 



TEMPERATURE CORRELATION 



Temperature fluctuations in the multiply-connected FRW cosmological models can be written as linear 
combinations (using £„; m in the last section) of independent components of temperature correlations in the 
simply-connected FRW cosmological models. Assuming that the perturbations are adiabatic and super-horizon 
scalar type and the initial fluctuations are random Gaussian, the two-point temperature correlation in a CH 
cosmological model is written as 



where 



^(„)^(„')\= ]T (a^a^YtMY^in*) 
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x ' i/(z/- 2 + l)Vol(S2) 

L v i m (r] ,n) = -Yi m (n)F l/ i(r] ), 

1 /" )o d$ t 

F v i{r) )=-$t{'n*) x vi{'no-'n*)+2j d v ~^ Xi,l< < rio ^ ri ^ 

5(sinh 2 77 — 3?7sinh?7 + 4cosh?7 — 4) 
t(jl)= (coshr?- I) 3 ' 

Here, V$(v) is the initial power spectrum, Vol(fi) denotes the volume of the CH space and $ v and &t(v) 
are the ^-component and the time evolution of the Newtonian curvature perturbation, respectively. 77* is the 
conformal time of the last scattering and r\ is the present conformal time. The diagonal elements (7 = V and 
m = m') give the approximate angular power spectrum Ci as 
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(2Z + l)Ci= J2 < l a '™| 2 ) 
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l^lml 2 !^! 2 . (8) 



It should be noted that the non-diagonal terms (either ( ^ (' or m ^ m') are not negligible for large angular 
(small 1) fluctuations in anisotropic models such as CH models. Therefore, constrains on the models by using 
only the angular power spectra Ci are not sufficient. On the other hand, this property can be considered as 
the "fingerprint" of the multiply-connectedness of the spatial geometry of the universe. 

Computation of highly-excited eigenmodes is still a difficult task since the number of the eigenmodes increases 
as k 3 and the number of the boundary elements increases as k 2 . In order to avoid these difficulties, we assume 
that £„/ TO 's are random Gaussian numbers and the variance is proportional to v~ 2 . Weyl's asymptotic formula 
gives the approximate values of highly-excited eigenmodes 

where N is an integer. We use pseudo-Gaussian random numbers £„; m that are derived from 14 eigenmodes on 
m003(-2,3) using the DBEM for k< 10, and random Gaussian numbers £,v(N)im whose variance is proportional 
to v- 2 for 10<fc<20. 

In FIGURE 1, ST/T t = [1(1+1)Ci/2tt] 1 ' 2 is (diamonds) plotted with the COBE data analyzed by Gorski (stars) 



[12] assuming that the initial power spectrum is the (extended) Harrison-Zeldovich spectrum V<s>{v) — Const. 
The 1-cr error bars are obtained by Monte-Carlo simulation with 10000 realizations. ST/Ti is almost constant 
in the limit o to 1 for the Harrison-Zeldovich spectrum. We see from these figures that ST/Ti is almost flat 
for 0.2 < f2 Q < 0.6. Suppression of the large angular power due to the long "wavelength" cutoff is quite mild 
compared with some fiat multiply-connected models since the bulk of the large angular power comes from the 
decay of curvature perturbations well after the last scattering time, which is known as the integral Sachs- Wolfe 
effect. Considering the cosmic variance, the suppression of the large angular power for the O = 0.1 model is 
still within the acceptable range. 

SUMMARY 

We numerically obtain 14 eigenmodes on a compact hyperbolic (CH) space m003(-2,3) with volume 0.98139 
using the direct boundary element method (DBEM) . The temperature fluctuations are written in terms of the 
expansion coefficients ^ u i m and eigenmodes on the universal covering space. For the 14 eigenmodes, £„/ m 's are 
numerically found to be pseudo-random Gaussian numbers with variance proportional to v~ 2 . 

The angular power spectra are computed using the 14 eigenmodes and an approximate method for eigen- 
modes with large k which is based on the assumption that the expansion coefficients are Gaussian random 
numbers. In contrast to multiply-connected flat models, the suppression of the large angular power is found 
to be so weak that the obtained powers are consistent with the COBE data for 0.1<f2 o <0.6. Assuming that 
the initial perturbations are adiabatic, constraints on CH models are not so severe as long as one uses only the 
angular power spectra which contain only isotropic information. 

However, one must also consider the anisotropic information of the temperature fluctuations. Contribution 
of the non-diagonal elements to the two-point temperature fluctuations is one of the key issues. Recently, Bond 
et al have obtained much severe constraints on the size of the topological identification scale for CH models 
using a method of images [13]. At the moment, the relation between their result and the author's result is not 
clear. Various methods for extracting the anisotropic information have been suggested such as a search for 
circles in the sky [6], or pattern formation [14]. The searches for the multiply-connectedness in the universe 
have just begun. 
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